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Abstract. We calculate the trispectrum of the primordial curvature perturbation 
generated by an epoch of slow-roll inflation in the early universe, and demonstrate that 
the non-gaussian signature imprinted at horizon crossing is unobservably small, of order 
t nl ^ r/50, where r < 1 is the tensor-to-scalar ratio. Therefore any primordial non- 
gaussianity observed in future microwave background experiments is likely to have been 
synthesized by gravitational effects on superhorizon scales. We discuss the application 
of Maldacena's consistency condition to the trispectrum. 
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1. Introduction 

In the inflationary scenario [H El El HI EJ El [7j , the primordial curvature perturbation, (, is 
generated via the vacuum fluctuations of one or more light scalar fields [HI El UOl [Hj . It is 
the only relevant perturbation in the simplest class of single-field models. More generally, 
however, ( is sensitive to the model under consideration. It may be accompanied 
by correlated or uncorrelated isocurvature perturbations [121 H3] which source its 
evolution, or it may be generated after the end of inflation by the decay of another 



It has recently become clear that interesting information about the particle physics 
that drove inflation is encoded in the non-gaussian features of the curvature perturbation 



[IS[T9l[20l[Hl[22l[2a[2l[2g 



HO] . Although observations of the Cosmic Microwave Background (CMB) anisotropy 
constrain the spectrum of perturbations to be close to scale-invariance with very nearly 
gaussian statistics [HI H21 021 HI] , some non-gaussianity is inevitable due to the universal 
coupling of all matter fields to gravity [151H61H7] . The question remains, however, as to 
the precise level of non-gaussianity that is generated. This is sensitive to details of the 
inflationary model, such as the number of fields which contribute to the energy density 
of the universe during inflation [HI H9J [501 EH E2], and also to possible non-canonical 
structures in the scalar field Lagrangian [531 EH ES, ESI E3 EE, EH] . 

If pure gaussian statistics are respected, any n-point correlator of the curvature 
perturbation ((■■■() can always be expressed in terms of the two-point correlator (C0- 
A breakdown of this rule for any n ^ 3 indicates a departure from gaussianity. Since 
some higher-order correlations may be more easily detectable [601 ED E2] than others, it 
is important to quantify the theoretical predictions for as many values of n as possible. 
To date, the majority of quantitative theoretical studies have focused on the three-point 
function, or equivalently the 'bispectrum' B^, 



Precise gaussianity corresponds to Bq = 0. In practise, the non-gaussianity in the three- 
point function generated in a particular model of inflation can be quantified in terms of 
a dimensionless quantity, /nl ; which is defined by [HI [631 [6lj|l 



where P^(ki) is the power spectrum of (, such that (C(ki)C(k 2 )) = (27r) 3 <5(^V kj)P^(A;i). 
Current observational bounds require |/nl| ^ 100 [65] and are consistent with /nl = 0, 
whereas a value of |/nl| ~ 3 should be detectable in the near future [53] ■ On the other 
hand, it is a robust prediction that the non-gaussian signal imprinted at horizon crossing 
during single-field inflation is |/nl| ~ 0.01, which is too small to be observable, [Tj5] and 
the situation is similar in models with more than one field [221 1251 HH [261 12H H3 150] . 
(For an alternative analysis, see Refs. [661 167].) 

| An alternative sign convention is used in Ref. [53] . 



field m Hang in]. 




(i) 



B^ = — /nl[-Pc(^i)-^c(^2) + cyclic permutations], 



(2) 
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In this paper, we make a comparable quantitative prediction for the primordial 
four-point function, or equivalently the 'trispectrum' [68, [69]: 



This can be quantified by defining a dimensionless 'non-linearity' parameter tnl such 
that [321 EH EE] 



where kj,- = kj + kj, and the condition = ensures that one-half of the 

permutations are equal to the other half, giving twelve distinct terms. 

At present there is only a weak experimental bound on the non-linearity parameter, 
which is roughly |tnl| < 10 8 [2H]. The WMAP satellite should strengthen this to 
I t nl| ^ 2 x 10 4 and the Planck satellite will be sensitive to a value of |tnl| ~ 560 [69J. 
Therefore, there is a pressing need to develop estimates of tnl in those models of inflation 
which we wish to confront with experimental data in the near future. An expression for 
r NL was presented recently by Alabidi & Lyth [28] (see also [31]). However, these authors 
neglected any contributions arising from quantum interactions at horizon crossing. The 
purpose of the present paper is to provide a more complete expression for the non- 
linearity parameter where such quantum effects are taken into account. We consider a 
general multi-field scenario in which the fields are minimally coupled to gravity with a 
target space metric 5 a p, and where the slow-roll approximation applies around the time 
of horizon crossing. 

The plan of this paper is as follows. In §2], we outline the 5N formalism for 
computing the trispectrum of the curvature perturbation on super-horizon scales. In §|3l 
we calculate the fourth-order interaction lagrangian for the field perturbations, including 
both scalar and vector contributions to the action for the first time. Given the form 
of this lagrangian, the four-point expectation value of the field fluctuations, and hence 
the momentum- dependence of the trispectrum, can be calculated by employing the 
Schwinger-Keldysh formalism [70j [7_H [72j [73] . This is achieved in §H]and represents one 
of the main results of the paper. We then derive an upper bound on the magnitude 
of the trispectrum in §5] and find that tnl < r/50, where r < 1 is the tensor- 
to-scalar ratio. Such a result implies that any non-gaussian signature generated at 
horizon crossing from the four-point correlator of the curvature perturbation will be 
unobservably small. We then consider the case of single-field inflation in §j6j where 
the curvature perturbation is automatically conserved on superhorizon scales. This 
conservation leads to a 'consistency' condition between the correlators of ( of the type 
discussed by Maldacena [19]. We conclude with a discussion in §[7J 

2. The 5N formalism for the trispectrum 

The most efficient method for calculating the four-point correlator of the curvature 
perturbation is to employ the 5N formalism introduced by Starobinsky [Til 175] . and 




(3) 
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generalized to non-linear perturbation theory by Lyth, Malik & Sasaki [TBjJi In this 
approach, it is assumed that the evolution of the universe in causally disconnected 
regions is like the evolution of separate locally unperturbed universes, where pressure and 
density can take different values [SUES ED!- This implies that C = 5iV[3llZ5lCSlE3ll211, 
where N(4>, p) is the number of e-folds of expansion between an initial spatially flat 
hypersurface on which the fields have values a and a final uniform-density slice on 
which the energy density is given by p. On large scales, ( is equivalent to the comoving 
curvature perturbation, 71, up to a sign convention. 

As demonstrated by Lyth & Rodriguez [33], if slow-roll inflation is valid on the 
initial flat slice, £ can be written as a power series in terms of the field perturbations at 
that time 

C = Njr + -N^S^ 13 + ■■■, (5) 

where iV jQ , = dN/d<f) a and a indexes the space of light fields <ft a . The perturbations in 
these fields are denoted by 5<j) a and are defined on the flat hypersurfaces. The leading- 
order relation for the power spectrum of the curvature perturbation is then given by 

P c = 5^N a N^, (6) 

where a subscript V implies that background quantities should be evaluated at horizon 
crossing and P* is related to the dimensionless power spectrum, V*, for a massless scalar 
field by 

V*{k) = k 3 P*(k)/2ir 2 = (Hj2n) 2 . (7) 

In principle, receives corrections from higher derivatives of N. We will assume that 
these corrections can be neglected when computing the leading non-gaussianity. This 
is reasonable since it is known experimentally that ( is dominated by a very nearly 
gaussian contribution. 

To calculate (CCCC) ? one forms the product of four copies of Eq. (jHJ) and takes an 
expectation value in the initial vacuum state [81j . The lowest-order contribution in 
derivatives of N is given by 

A(C(k 1 )C(k 2 )C(k 3 )C(k 4 )) = N^N.Nsidr^i^i^^^k^^)), (8) 

where the correlator of field perturbations on the right-hand side is evaluated a few 
e-foldings after horizon crossing, such that the fields have had sufficient time to become 
classical but have not evolved significantly. The 50-correlator in (jHJ) can be separated 
into two pieces. The first is an irreducible connected part in the sense of Feynman 
diagrams (see Figure Ufa)), which we denote by ((50 a 50 /3 (50 7 ^0' 5 } c , and which is absent 
when the fields perturbations are precisely gaussian. The second contribution is a 
reducible part and is given by the sum over all ways of combining the four fields 
into pairs, with each pair yielding a copy of the spectrum. This reducible part is a 
disconnected contribution, and is always present, even when the 5(f> a are gaussian. It 
therefore contains no more information concerning the primordial non-gaussianity than 

§ See also Langlois & Vernizzi |77j for an alternative non-linear generalization. 
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Figure 1. Connected (a) and disconnected (b) contributions to the four-point function 
of the {5(f)} . The disconnected contribution factorizes into a product of two copies of the 
power spectrum, and contains no information concerning primordial non-gaussianitics 
beyond those encoded in loop corrections, which are expected to be unobservably small 
(821. 



is already present in the power spectrum, and in the remainder of this paper we consider 
only connected contributions to the trispectrum. In principle, the contribution to (CCCC) 
in Eq. (|Sj) is accompanied by terms which are higher-order in the derivatives of N, of 
which an example was computed by Sasaki, Valiviita & Wands [83]. We will briefly 
address these contributions in and will return to them in a forthcoming publication. 

In obtaining a prediction for the non-linearity parameter t^l from Eq. (jSJ), it proves 
very useful to decompose the trispectrum into a number of copies of the dimensionless 
power spectrum, 7-*. For a target space metric 5 a p, one may define the connected part 
of the trispectrum in terms of a form-factor Ai^hi, k2, k3, k^, which parametrizes the 
momentum-dependence of (5(j)5(j)5(j)5(j)). More specifically, combining Eqs. (jSJ), © and 
(jEJ) leads to a definition for M.± given by 

T H\T = ^ 3 E SWM&u k 2 , ks, k 4 ), (9) 

1 L% i perms 

where i G {1,2,3,4} and the sum over permutations includes all rearrangements of 
the indices {a, ft, 7, 8} which simultaneously permute {ki, k 2 , k 3 , k 4 }. Note that the 
form-factor may depend on the relative orientation of the kj and not merely on their 
magnitude. Note also that the ki should have approximately equal magnitude in order 
for the epoch of horizon crossing to be more or less unique. 

The momentum-dependence of tnl may then be deduced by comparing Eqs. (jSJ) 
and flSj) with Eq. This is facilitated by rewriting Eq. in the form 

T c , c = 47r6 r NL (^N a N^V! £ E k " k % k ^ + ( 10 ) 

llie{l,2,3,4} K i i<j m^i,j 

where we have substituted for the power spectrum ([6]). It then follows from Eqs. (jSJ), 
(EJ) and (TTUTl that the contribution to tnl that is generated at horizon crossing is given 
by 

a E pe rm S -^4(kl,k 2 ,k 3 ,k 4 ) 

TNL (5-m, a N,)E t<J E m ^kfk^ + kjiy [ j 

where the sum is over all permutations of the {ki, • • • , k^. 
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Finally, the non-linearity parameter can be related to the 'tensor-to-scalar' ratio, 
r, which is defined in terms of the corresponding power spectra such that 

where Eq. ([6]) has been employed. Hence, it follows that 

r Y M A 

A _ Sperms'' 4 /i o\ 

NL ~ 8V V k 3 k 3 (k~ 3 + k~ 3 Y { } 

° l^i<j 2^m^i,j K i \ K im + K jm) 

The magnitude of this contribution can therefore be derived immediately once the form- 
factor A^4(ki, k2, IC3, ki) has been calculated [25J. This calculation will be the focus of 
the following two sections. 



3. The fourth-order action for the perturbations 

The appropriate formalism for computing expectation values such as (5(j) a 5(fiP Sep 5(j) s ) 
is the so-called 'closed time path' integral, which was developed by Schwinger [70J and 
others [7T| |gl]. and subsequently extended to cosmology [723 ESI ESI ESI ESI ESI E7J • This 
integral expresses the four-point correlator as 

(54> a (t', kJS^it', k 2 ) W, k 3 )5(P\t', k 4 )> = 

j [d«50 + d50_] ^(t / ,k 1 )^(t',k 2 )50X(t',k 3 )^ + (t , ,k 4 ) expiH/, (14) 

where W is a weight involving the action for the field perturbations 5<f) a and takes the 
form 

W[5(p + ,5(f).}= j dtd 3 x L{5<f) + ) - j dtd 3 xL(50_), (15) 

where L(d~<j)) is the lagrangian for the {6(f)}. In Eqs. (|T4l) - ffT5l) . we have adopted a 
common time of observation t', and is any fiducial time satisfying t' < t*. The 
functional integral is over all 'forward-going' fields 5<j) + [H] which begin in the vacuum 
state at t — > — oo together with all 'backward-going' fields which obey the same 
boundary condition and coincide with the S(p + at t*. This differs from the conventional 
Feynman expression in which only the forward-going field is present [88| [89] . This 
difference arises because the Feynman formula computes a scattering amplitude from 
the past to the future, whereas the Schwinger formula computes an expectation value. 

The lagrangian can be naturally expressed as a sum of terms involving a definite 
power of the perturbation S(p, 

oo 

L{5cj ) ) = Y,L k m (16) 

k=2 

where contains k such powers. Hence, the gaussian free theory is completely specified 
by L2, but the Lk terms with k ^ 3 correspond to interactions. In the perturbative 
regime these can be assumed to be small corrections to the free theory, and generate 
the non-gaussianities we wish to calculate. In order to compute the n-point expectation 
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value, it is necessary to know the functional forms for all with k ^ n. The L 2 
and L 3 terms were obtained in Ref. [22] and are sufficient to calculate the three-point 
expectation value and therefore the bispectrum Bfy. The term L4 is necessary to 
determine the four-point expectation value, and hence the corresponding trispectrum 
Ts^. The scalar part of L 4 has already been derived by Sloth jH2]. In the following 
section we extend this result by incorporating the effect of the vector part of L 4 . 



3.1. The Lagrangian for field perturbations 

We consider the class of inflationary models where Einstein gravity is minimally coupled 
to a set of scalar fields <p a : 

S = Jdtd 3 x^g~ Qi? - i V> Q V a Q - V^j , (17) 

where g = detg, R is the Ricci scalar of spacetime and V is the interaction potential 
for the fields. It is assumed that the potential can drive a phase of slow-roll inflation, 
but is otherwise left arbitrary. 

In order to work with the 5N formalism, we require the perturbations 5<p to 
be defined on 'flat' spatial hypersurfaces on which the three-metric takes the form 
hij = a 2 (t)[5ij+exp(jij)}, where d l ^ij = 0. These slices are indeed flat in the conventional 
sense if the tensor mode 7^ is absent, which for simplicity we assume to be the casqjj. 
It proves convenient to employ the ADM form for the metric, 

ds 2 = -N 2 dt 2 + hi^dx* + N { dt)(dx j + N j dt), (18) 

where N(t,x) and N l (t,x) denote the lapse function and shift vector, respectively. The 
action (jPTl) can be expressed in terms of these degrees of freedom such thal0 

S = -^Jdtd 3 x NVh (V> a Vi0 a + 2V) + ^ J dt d 3 x ^ {E^E^ - E 2 + 7r Q vr Q ) , (20) 

where Eij = \h%j — V(jiVj) is proportional to the extrinsic curvature of the spatial slices 
and 7r a = <p a — N^Vj4> a is proportional to the momentum conjugate to <f) a . 

The degrees of freedom iV and iV* may be integrated out by first extremizing the 
action with respect to these fields, solving the corresponding algebraic equations of 

|| On the basis of the calculations reported by Maldacena in Ref. [TH], any expectation value involving 
a tensor mode will be smaller than an expectation value involving only scalars. Since the scalar non- 
gaussianity is already hard to detect (and, despite a strong theoretical prejudice in its favour, the tensor 
contribution has not yet even been confirmed to exist), it seems reasonable to restrict attention to the 
scalar sector. In that case, 7^ does not enter in the tree-level graphs and can be dropped without loss 
of generality. 

In what follows we employ a useful summation convention for the spatial indices, which are labelled 
{i, j, ■ ■ •}. When these appear as a raised/lowered pair, it is to be understood that they are contracted 
with the 3-metric, hij. On the other hand, if a pair appear with both indices lowered, they are 
contracted with the flat Euclidean metric 6ij. In other words, 

a l bj — ^""^ h^ajbj and ajbj = ^""^ ajbj. (19) 
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motion, and then substituting the solutions back into the action. This results in an 
action that can be expressed entirely in terms of propagating fields. Although the 
solutions to the constraint equations are in general non-local, which is signalled (for 
example) by the presence of formal operators such as the inverse Laplacian d~ 2 , this will 
not be problematic provided we remain interested only in momentum-space expectation 
values. On the other hand, the continuation of such expectation values back to real 
space may not necessarily be a straightforward procedure. 



J2 a m, $=Y1 & = Yl P<»3> ( 24 ) 



3.1.1. The constraints The constraint equation for the lapse function is given by 

V'fV,^ + 2V + {EVEy - E 2 + n a n a ) = (21) 
and the corresponding expression for the shift vector takes the form 

(^m - m) = ^oK. (22) 

Eqs. (121]) and (1221 can be solved perturbatively by writing 

iV = 1 + a, Nj = + Pj, (23) 

where (3j is diver genceless, i.e. Vj/3 J = 0. It is then assumed that the variables a, d 
and Pj can be expanded as a power series in the field perturbations: 

oo oo 

a = 

m=l m=l m=l 

where a subscript m denotes the number of powers of 5<p that are present in each term. 
As shown in Ref. [Mj, although in principle we would need to include terms in the sum 
up to m = 4 in order to obtain L 4 , in practise the terms with m = 3 and m = 4 cancel 
out of the final answer. 

Let us first consider the N l constraint. We find from the 0(50) term that 

ai = ^rS(f> a and ^ = 0, (25) 

whereas a2 and faj are obtained from the 0(5<ft 2 ) term: 

2HdjCt2 — -^d 2 ^ = iHaidjai + 54> a dj5<f) a — ai<fi a dj5<f) a 

+ \d i a l (5 ij d 2 $ 1 -d i d j $ 1 ). (26) 

In principle, the O(<50 3 ) term could be used to find a 3 and j3 3 j, but these contributions 
are not needed. 

With regard to the lapse function, the 0(1) term in the constraint equation (|2~T!) is 
the Friedmann equation: 

6H 2 = 2V + <p a <p a . (27) 

The 0(50) term yields $i, 

i^«9 2 ^i = -2V a 54> a - 2<p a 5<p a + 2aU-6H 2 + a 0A (28) 
a 2 
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and the term $ 2 follows from the 0{5<p 2 ) equation: 

%dH 2 = -Id^rd^ - v af3 54> a S4>P 

a z a z 

+ \d 2 ^d 2 ^ - - 5<j> a 5<j> a + 2^diMl 



a ar a* 

V a 2 



+ 2ai ( — d 2 ^ + 2<p a 6<p a ) - (3a 2 - 2a 2 )(-6H 2 + 0> a ). (29) 



Likewise, the 0(5<f) 3 ) term can be used to obtain $ 3 . Note that since !? 2 / and it 
remains finite as <p a — > 0, some metric back- reaction from the coupling to the scalar 
fields will still be present even in the extreme de Sitter limit. 

The second-order perturbation theory of multiple scalar fields has previously been 
studied by Malik in the large-scale limit [90]. Compared to the lower-order perturbation 
theory considered in Refs. [19], [531 122], the most distinctive new feature of the above 
analysis is the combination of a 2 and (5 2 j appearing in Eq. (|26l) . These terms can be 
separated using the divergenceless condition on (3j together with the reasoning that 
is employed in the proof of the Helmholtz Theorem [9TL [92] . For this purpose, it is 
convenient to recall the elementary identity 

V 2 W = grad div W - curl curl W, (30) 

for any vector field W. A given vector F may be decomposed into its incompressible 
and irrotational pieces by writing F = grad cp + curl A and using Eq. fl30|) to identify 



1 1 
A = -curl— F and <p = div — F. (31) 

After applying this construction to Eq. (|26|) . it follows that a 2 can be expressed in 
the form 

oii = 4 + T77 d ~ 2 f S + \d 2 ai d 2 A - -d^a^d^ , (32) 
Z An \ a a J 

where 

£ = djd^djd^ + 5<p a d 2 5(j) a . (33) 
Likewise, it can be shown that to leading-order, (3 2 j satisfied 

^/3 2j ~ «T 4 («9 m «V0 Q a m 50 a + d 3 5^ a d 2 5<p a - d 2 5rd j 5ct )a - «9 m 50 Q «9 m <V0 Q ) (35) 
where '~' denotes equality only at leading-order in the slow- roll expansion. 
+ The exact expression is 

r^Aj = <9~ 4 ( \d 2 aidjd 2 i9i - \d m dja\d m d' 2 'dx + -\d m aid m djd 2 $i - \djaid i '&i 
2a 2 \a 2 a 2 a 2 a 2 

odmdjdiaididm'&i + -^d 2 d i aid l djdi ^djdiaidid 2 t)i + —zd m diat\d m didj'&i 

a 2 a 2 a 2 a 2 

+ d m ,d j &j> a d m 64> a - d 2 S<j> a d j Sct) a + d j 84> a d 2 5(t> a - fl^ a flU0 3 -ty a ) ■ (34) 

It is important that foj cannot be written as a total gradient. Similar terms were identified in Refs. 
[231 177) . and in Ref. [77) it was shown that such quantities decay like a~ 3 in an expanding universe. 
This is in agreement with non-perturbative arguments which predict the decay of shear during inflation 
[93] . and is essential to maintain the relationship between TZ and C on large scales. We thank F. Vernizzi 
for emphasizing this point to us. 
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3.1.2. The fourth-order action After use of the constraint equations §ZT§ and 
outlined above, the fourth-order action is found to be given by 



as 



+ (ala 2 --a 2 2 ){-6H 2 + 



Ha) + ^ 



- 2alV 



+ Oil 



V, 



- - i d i d j 2 did j <& 1 + -d 2 $ 2 d 2 A - -difajdidfl! 
a a Qj 



(36) 



dj'&idji 

a- ' ' " a' 

When appropriately truncated in powers of slow-roll, the scalar sector of this action 
is equivalent to that presented in Appendix A of Ref. [82]. However no truncation in 
slow-roll has been made in Eq. (136]) . which is exact. 

The leading-order, slow-roll terms in Eq. ( 1361) arise from a 2 , d 2l and (3 2 j [via 
Eqs. (J32J), flU and ([35}]. These all contain terms of 0(e°), where e = -H/H 2 
parametrizes the degree of departure from a pure de Sitter expansion. Extracting these 
contributions from Eq. f[36l) then leads us to the fourth-order action for the perturbations 
that we have been seeking: 

1 „ „n „ _ _ 3 



S 



dtd'x ( -—(3 2j d 2 (3 2j + m? 2 £ + -a J 9~'S^ 2 S - a^ a (5 2j d^ c 



(37) 



The form of the lagrangian L4 then follows immediately from Eq. (1571) . 



3.2. Loop corrections 

The presence of interaction terms in L4 which are 'unprotected,' in the sense that 
they are not suppressed by a positive power of the slow-roll parameters, might give 
rise to concerns that the two-point function (<50 a 50 /3 ) also receives large, unprotected 
loop corrections. If this were the case, it would seriously impair our ability to make 
predictions about the early universe from models of inflation. Fortunately, however, the 
calculations reported in Ref. [82] indicate that these unprotected loop corrections occur 
at about 1% of the tree-level, which is not large enough to harm our predictivity. On 
the other hand, it should be remarked that it has yet to be verified that such corrections 
are never large at any order in the interaction of 8(f). 

Given the interaction lagrangian (1371) . it only remains to apply Eqs. (TT4"I) - (115D in 
order to obtain the four-point expectation value of the field fluctuations and therefore 
the momentum- dependence of the form-factor M.4 and the non-linearity parameter tnl- 
This will be focus of the following section. 



The inflationary trispectrum 



11 



4. The trispectrum form-factor M. 4 

4-1. The momentum- dependence of M.4 

It is convenient to calculate the four-point expectation value in three separate stages 
by considering the blocks of terms in Eq. ( 1371) which consist of either two, one or zero 
copies of the pure vector (3 2 j- 

• Vector x vector. This is the term containing (3 2 jd 2 (3 2 j. It can be calculated 
following the general principles outlined in £J3] and described in more detail in 
Refs. P21IS51I531I221IB51E21CZ21CZ31E5. Before proceeding with the calculation, 
however, it will be helpful to introduce some new notation. In particular, we define 
a vector which is a function of two momenta k« and k,, such that 

Zjj = <7ykj — cjjikj (no summation), (38) 

where is the scalar combination 



a 



'j 



ki ■ kj + k). (39) 



The vector x vector term then provides a contribution to the four-point expectation 
value (Scf) a (rj', k 1 )50 /3 (?7 / , k 2 )<5</> 7 (7/, k 3 )<50 5 (7/, k 4 )) which is given by 



A ±±2 2 



e-W^^^Zu ■ Z 34 



/n* 
dr] H(r)) 2 (l - ih 2 T))(l - ik A rf)e lktr) + permutations + c.c, (40) 
-00 

where 'c.c' denotes the complex conjugate of the preceding term, the index i 
ranges over {1, 2, 3, 4}, the conformal time equivalents of (£', £*) are given by (?/, 77*) 
and the total scalar momentum is defined by k t = k\ + k 2 + ^3 + k A . The sum 
includes all permutations of the indices {a, j3, 7, 5} which simultaneously permute 
the momentum labels {1, 2, 3, 4}. 

The time of observation, 77', should be chosen to be a few e-foldings after the modes 
kj have crossed the horizon, so that it can be safely assumed that the fluctuations 
5(f) have become classical but have not undergone significant evolution. At the 
level of accuracy to which we are working, this will be the case if the condition 
\k t f] f \ ~ 0(e) 1 is satisfied. When this condition holds, we may further specify 
if (7/) = H(i]) = to a good approximation, since the integral over 77 will receive 
only negligible contributions from times long before the epoch of horizon crossing. 
Likewise, extending the upper limit 77* in Eq. (j4"0]) to the infinite future (77* — > - ) 
introduces no significant error. This implies that the integral ( I4"0l) can then be 
evaluated analytically and we find that the contribution of this piece to the form- 
factor Ml 4 is given by 

AA* 4 = -2-gM_^Z 12 .Z 34 , (41) 

fc 12 fc 34 Kf 
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where we have defined 



K) 



ik f 



drj (1 — ikiT])(l — ikjT])e 



kj ~\~ kj ^kjkj 
1 + — — - H -y~ 

k t k t 



(42) 



Since Eq. ( 1421) is purely real, the addition of its complex conjugate has been 
absorbed into an overall factor of 2. 

Scalar x vector. The same principles apply to the scalar x vector sector, which 
comprises the term a5(p a P2jdjS4> a . Its contribution to the four-point expectation 
value is given by 



ft 34 



krZ 



31 



dr] (1 — i/c2^)(l — ik^rj)e v + permutations + c.c. 



(43) 



Since the integral is equivalent to W24, we may immediately deduce the 
corresponding contribution of this term to the form-factor is 
k\k\ W 2 



AA4 4 = -4- 



'24 
^34 "4 



k2 ■ Z34. 



(44) 



Scalar x scalar. The final piece is comprised of + (3/4)a 3 <9 _2 E<9 _2 E and is 
independent of f3f. A loop integral involving this term was calculated in Ref. |82j . 
After truncating $ 2 to leading order, using Eq. (f2"9|) . the scalar piece can be rewritten 



'scalar 



dt dx 



AH 



-<T 2 E 



; di5(j) a di5(f) a - 



- -a 3 cr 2 E(T 2 £ 
4 



(45) 



The first of these terms makes a contribution of the form 

1 A; 2 



x 



dr? [ki • k 2 (l - i^)(l - ik 2V ){l - ik 4V ) - k\k\tf(\ - ^477)] e ifct " 
+ permutations + c.c. 



(46) 



It proves useful when evaluating this integral to define a i/iree-momentum 
generalization of Wij, 

"0 



Iran 



ik t ri 



drj (1 - ik e r))(l - ik m r))(l - ik n r))e 



kg + k m + k n 1(k(k m + k(k n + k m k n ) Qk(k m k ri 
1 , 1 , o ~ r 



Hence, the contribution of Eq. (1461) to the form-factor is 



1 k 2 

AM 4 = -^Tf-0- 34 0- 34 
z fc 34 



: Wl2A H 



2 + 6 



Finally, the second term, — (3/4)a 3 <9 -2 £<9 -2 £, yields 

ff6 3 1.2 u2 



(47) 



(48) 



, , „.., kiki 
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x 



/drj (1 — ik 2 rj)(l — ik^rf)e lkt ' n + permutations + c.c. 
-oo 



and provides the contribution 

3 A?' 



AM 



4 - -^rjr°^ a ^-L 



W; 



Hi 
t 



13 

(49) 

(50) 



to the form-factor. 



It follows, therefore, that the complete momentum-dependence of the trispectrum 
can be determined by assembling Eqs. (l4Tj) . (jUJ), (PfB]) and (!50|) . Specifically, we find 
that the form-factor is given by 



M 4 



, kjkj W 24 

' k 2 k 2 k< 



z 12 ■ z 



31 



^34 



2k 2 • Z 



31 



^C"l2C"34 



1 A; 2 
-— (T 3 4 



2 fc 



31 



6- 



— ' fc? V~ 1 "WJ' (51) 

and Eq. (IBTj) represents the main result of the paper. The form-factor has dimensions of 
three powers of the momentum, which is the same behaviour as the corresponding form- 
factor which determines the bispectrum. However, the momentum-dependence of 
M.4 is considerably more complicated than that of M.3 and, in particular, it depends 
on the products kj • kj and not merely the magnitudes fcj. 



4-2. The limit of equal momentum 

If two of the momenta, say k x and k 2 , have equal magnitude k\ = k 2 = p, and are also 
equal and opposite, kx + k 2 = 0, it follows from momentum conservation that the other 
two momenta must also be equal and opposite, albeit with a possibly different common 
magnitude, i.e., fc 3 = k A = q. On such configurations, M.4 appears to be pathological, 
since the factors (ki + k 2 ) _1 and (k 3 + k 4 )~ 2 naively diverge. If this were the case it 
would be disastrous, since the non-gaussianity associated with the four-point function 
is presumably bounded from above by the WMAP data. 

Despite appearances, the form-factor M.^ does not diverge on these degenerate 
quadrilaterals. This can be verified by introducing a simple specification for the kj. We 
define a matrix of angles 6^ G [0, 7r] which satisfy 

k ■ k 

cos% = 13 for i ^ j. (52) 

This matrix obeys the obvious symmetry constraint % = 9ji, and therefore has 6 
components, which cannot all be chosen independently. To fully specify the {kj}, it is 
necessary to include two additional angles which determine the absolute orientation of 
one of the kj, which we take to be k 4 . However, these angles do not contribute to .M4 
and we need not be explicit about their precise assignment. There is another angular 
degree of freedom corresponding to rigid rotations which leave k 4 fixed, which does not 
change the intrinsic geometry of the {kj}. In total, therefore, the magnitudes of the 
four momenta, the two angles which determine the orientation of k 4 , the angular degree 
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of freedom corresponding to rotations which leave k 4 fixed, and the five independent fly- 
comprise the twelve degrees of freedom which are required to completely specify four 
vectors in R 3 . 

Momentum conservation supplies three constraints which eliminate three of these 
degrees of freedom. These can be chosen to be three of the angles fl^-. It is convenient 
to work with the angles {fl r = fl r4 }, in terms of which the others can be written as 

cos 9 13 = cos fl 2 , (53) 



cos fl 23 

COS 9\2 



COS fli , 
2 



Q_ 

p' 2 



;i + cosfl 3 ) - 1. 



(54) 
(55) 



We wish to re-express the form-factor (|5"T|) in terms of p, q and 9 r . This requires 
explicit expressions for TV24 and W124, which can be shown to be 

6(p + q) 2 + 2pq 



24 



4(p + g) 3 



(56) 



8(p + g) 3 + 4(p + g) 2 (2p + g) + 4p(p + g)(2g + p) + 6p 2 g 
^ 124 " WTqT • (5?) 

We also have the additional relation o~ 12 = cr 34 = g 2 (l + cosfl 3 ) and the quantities Z 12 -Z 34 
and k 2 • Z 34 can be reduced to 

Z 12 ■ Z 34 = 2pq 5 (l + cosfl 3 ) 2 (cosfl 2 - cosfli) (58) 
k 2 • Z 34 = ~pq 3 (l + cosfl 3 )(cosfl 2 - cosfli), (59) 

respectively. 

It is now relatively straightforward to rewrite expression ( !5T|) given Eqs. (!56l) - (l59l) 
and we find after some algebra that 



Mi 



W- 



21 



2 

p q 



4(1 + cosfl 3 )p + q 
q 2 p 2 



3q 

p(cosfl 2 — cosfli) — + cosfl 3 ) 



8(p + q) 



W 



121 



P 



1 + cosfl 3 



1 + 



p 2 (2p + 5q) 



(60) 

7r. However, 



4(p + q) 3 

The naive divergence associated with the limit \ti — > — k 2 reappears as fl 3 
Eq. fl60l) should not be considered in isolation, but only after symmetrization over the 
labels {1,2,3,4}. Indeed, after symmetrization over the exchange 1 ^ 2 we find 



2^2 



MAI 



3p q 



-W< 



p 2 q 2 



ijp + q)" 24 4(p + g) 



w 



121 



%(l + cosfl 3 )-l 



p 2 (2p + 5q) 



.(61) 



4(p + g) 3 

and it follows that Eq. ( l6Ti) possesses a finite limit as fl 3 — > 7r which is not spoilt by 
further symmetrization. Hence, as anticipated, Eq. flBTT) demonstrates that .M 4 does 
not diverge on degenerate quadrilaterals. 
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5. How large can the trispectrum be? 

We now proceed to estimate the magnitude of the non-linearity parameter |Atnl| that 
arises from Eq. (15"TT) . This is important, since the non-gaussianity associated with the 
trispectrum may be observable in the near future if |Atnl| > 560 [69] . 

The form-factor M.^ is formally valid only when ki ~ kj for all i and j. We will 
show in §|6]that, at least in the single-field case, .M 4 must become of order 0(e) in the 
limit where one kj becomes much smaller than the others. In this limit one cannot trust 
( l5Tj) to give the dominant contribution to .M 4 , but it is consistent with the explicit 
momentum dependence in ( 15~TT) . which approaches zero in the formal limit where any 
ki — > 0. This provides evidence that the largest contribution to r NL from the connected 
four-point function of the {5<p a } arises on equilateral configurations, although a more 
complete treatment of this question would be desirable. In this section, we assume that 
|Atnl| is maximized on equilateral configurations of the {kj}. 

The relevant expressions for the equilateral configuration follow after we have 
specified p = q = k in Eqs. ( I53|) - (l59|) . The equilateral condition is not independent 
of momentum conservation. In practise this means that the angles 9 r cannot be chosen 
freely, but must instead obey the constraint 

^cos0 r = -l. (62) 

r 

We may choose only two of the 9 r arbitrarily. These two angles suffice to parametrize 
the family of distinct equilateral configurations. If k 4 is chosen to point along z, and we 
take (162]) to determine #3 in terms of {#1, 9 2 }, then this family can be written 

ki = (sin 6>i cos <px, sin 9\ sin cf>i, cos #1) 
k 2 = (sin 6*2 cos(0i + a), sin 9 2 sin(0! + a), cos 9 2 ) 
k 3 = (sin 6*3 cos(0i + (3), sin 6*3 sin(0i + ft), cos 6*3) 
k 4 = (0,0,1) 

where cos^ = — 1 — cos#i — cos 9 2 , and the angles a and (3 obey 

tana = ±A ^ an d tan/5 = T~ -, (64) 

73 72 

given that 7 m is a combination of the 9. 



(63) 



cos^ m -n r ^ m cos^ r 
7m = i— Q • 65 

The extra degree of freedom, 0i, corresponds to rigid rotations of the {ki, k 2 , k 3 } which 
fix k 4 cx z. One can verify explicitly that Eqs. fl53|) - fl55|) and Eq. fl62l) are obeyed by 
( l63l) . irrespective of the value of 4>x- In the equilateral limit, W 2 ± and Wvii reduce to 

1 S 

W« - y (66) 
We m n (67) 
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respectively, and this implies that Eq. (I6"0l simplifies to 



M 4 



T 



13 



1 



161 + cos 9? 



COS 02 — COS 6*i 



COS $3 



71 7 

COS C7q 

128 3 128 



(6* 



We may then substitute this result into Eq. (fT3l) . Summing over the possible 
permutations of labels and using Eq. ( 1621) to eliminate the sum ^ r cos# r , we deduce 
that 

23\/2r 

y (1 + cos 6 r 



An 



NL 



32 



\-3/2 



(69) 



The magnitude of tnl which is generated depends on the angles {9 r }. However, this 
dependence is in some sense artificial since the totally permuted form factor, Sperms -^ 4 > 
contains no angular dependence in the equilateral limit. The dependence on the {6 r } 
which appears in Eq. (1691) arises solely from the definition of tnl which we adopted in 



Eq. 



|Atnl| is maximized on the configuration cos# r 
23 r 

|Ar - l< wr^ - 023r ~5o- 



— 1/3, which gives 



(70) 



Since the current experimental upper bound on the tensor-to-scalar ratio is r < 0.55 (at 
95% confidence) [53], Eq. (170|) immediately implies that the primordial non-gaussianity 
in the trispectrum that is generated at horizon crossing will be unobservably small. 
Consequently, as in the case of the bispectrum, if a non-trivial primordial trispectrum 
is observed in a future CMB experiment, its origin is unlikely to arise from quantum- 
mechanical interference around the time of horizon crossing. 

We observe that ( 1701) is not zero. Taken together with the explicit family of 
equilateral configurations (163]) . which are all allowed by momentum conservation, this 
removes any worry that the complicated expression obtained for Ai^ in Eq. (13Tj) might 
be constrained to zero by purely kinematical considerations. 

Thus far, our discussion has considered the multi-field inflationary scenario. A 
number of simplifications arise when the analysis is restricted to single-field inflation 
and we proceed to discuss these in the following section. 



6. Single-field Inflation and the Maldacena consistency relation 

In single field inflation, ( is conserved to all orders on superhorizon scales [TBI [76] . This 
has the immediate consequence that no evolution in t^l is possible after horizon exit, 
and Eq. ( 1701) represents the largest possible non-gaussian signal which can be visible in 
the trispectrum. We conclude that it is too small ever to be detected. 

A second special feature of single-field inflation is the existence of a consistency 
relation between the correlators (C(ki) ■ ■ -C(k n )). In the 'squeezed' limit, where k\ — > 0, 
the C(ki) mode crosses the horizon at a much earlier epoch than the remaining modes 
{k2, . . . ,k n }. By the time these modes eventually exit the horizon, the gravitational 
background has been deformed due to the presence of the C(ki) mode. As first pointed 
out by Maldacena [19] . the ki mode and the {k2, ■ ■ • , k„} modes no longer interfere when 
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they cross the horizon, which implies that the only correlation between them is the one 
imposed by this gravitational relationship. In the case of the three-point function, this 
yields the relation 

(C(k 1 )C(k 2 )C(k 3 )) - -K - l)Pc(M(C(k 2 )C(k 3 )> (as h -> 0), (71) 

where n s is the spectral index of the scalar perturbation power spectrum. It was later 
emphasized by Creminelli & Zaldarriaga [94J that the limit (1711) is purely kinematical 
and applies for any metric theory of gravity. This implies that it can be calculated on 
the basis of gravitational physics alone |27j . 

Maldacena's arguments apply for any correlation function of (, including the four- 
point function. In this case, the consistency relationship is given by 

(C(k 1 )C(k 2 )C(k 3 )C(k 4 )) - -P(fc 1 )-^(C(k 2 )C(k 3 )C(k 4 )) (as h - 0). (72) 

To gain further insight into the nature of Eq. (175]) . it is instructive to count powers 
of slow-roll. In general, the derivatives of the number of e-foldings, N, are controlled by 
the slow-roll parameter e = —if /if 2 . If e is almost constant over the range of e-folds 
under consideration, it follows that N ~ — e -1 In if. This implies that each derivative of 
iV with respect to the scalar field generates an extra power of e 1//2 , i.e., ~ 0(e -1 / 2 ), 

0(e (n - 2)/2 ), etc. By counting powers of slow- roll in Eq. (172]) . 
it may be verified that this behaviour produces a four-point function of order 0{e^ 1 ). 
This can be compared to the horizon-crossing contribution, Eq. (jSJ), which is 0(e~ 2 ). 

Eq. ( 172]) implies that (|5T]) becomes at most of 0(e) when any of the ki — > 0, at least 
in the single field case. Therefore, the leading slow-roll term in the four-point function 
of ( will be given by possible subdominant corrections to (I51T) which do not vanish in 
the squeezed limit, together with the terms which contribute to the four-point function 
at next-order in the slow-roll parameter e, i.e., the terms of order 0(e _1 ). 

The first such term is 

(2tt) 3 5(^ k^-N^N^N^ [P{k 2 )B{k 12 , k 3 , fc 4 ) + 23 permutations] , (73) 

i 

where the sum contains twenty-four terms that are obtained by rearrangements of the 
momenta, and B is the <50-bispectrum. The momentum conservation condition implies 
that one-half of the permuations are equal to the other half, giving twelve distinct terms. 
The second relevant term is 

(2n) 3 5(J2 **)\ N ,W N ,<M> N ,* N ,+ [P{hs)P{h)P{h) + 23 permutations] . (74) 

i 

At the tree-level, only this term was included in the estimate for tnl given by Alabidi & 
Lyth [28] and Lyth [3l] , although these authors also included a one-loop diagram whose 
contribution we assume to be negligible in comparison with the tree-levelo 

* Owing to a typographical error, the tnl which follows from (|74[) appeared incorrectly in the journal 
version of Ref. The correct expression was given in both vl and v2 of Ref. [5S], available from 

the arXiv. [D. Lyth, private communication (2006).] 
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The final contribution at 0(e 1 ) is 
(2n) 3 5(J2 ^)\n,H4> N ,4> N ,4> N ^ [P (h) P (k 3 ) P (k 4 ) + 23 permutations] . (75) 

i 

An expression equivalent to this was given in Ref. [83] in the context of the curvaton 
scenario. 

In principle, therefore, the next-order contribution to the four-point function can 
be obtained by summing Eqs. ( J73l - (|75l) . However, when performing such a sum, self- 
consistency would require that the next-order versions of the two- and three-point 
functions of the field fluctuations 5(f) and the derivatives of N be employed. To date, the 
required expressions for these quantities have yet to be calculated since presumably they 
are unobservably small. Consequently, we should not expect the e and r\ dependence 
of Eqs. (I73l)- (l75l) to match the right-hand side of Eq. ( 1T2l . where n = V^/V is 0(e) 
in the slow-roll expansion. On the other hand, the next-order corrections to the scale 
factor, a(t), the wavefunctions of the 5<ft and the third-order action are not expected to 
contain any intrinsically second- or higher-order slow-roll parameters, which involve the 
third- or higher-derivatives of the potential. Therefore, summing Eqs. (173|) - (!75l) with 
the leading-order expressions for N, the scale factor a, and the two- and three-point 
functions of the field flucations 5(j) should match the coefficient of in Eq. f!72|) . We 
have verified that this is the case. 

7. Conclusions 

In this paper, we have derived the trispectrum of the primordial curvature perturbation 
imprinted at horizon crossing during inflation. In particular, we have estimated the 
expected magnitude of the non-linearity parameter, tnl, which is generated by this 
process. We find that it is bounded by about r/50, and therefore is too small to be 
observable by the Planck satellite or other future CMB experiments. In single-field 
models, ( is conserved to all orders on superhorizon scales. Therefore, no subsequent 
generation of non-gaussianity is possible, and the primordial signal will be too small 
ever to be detected. 

One unexpected feature of the four-point action is that it contains no powers of 
slow-roll parameters. In an earlier publication [22], two of us (DS and JEL) suggested 
that the nth order action would contain essentially \n/2\ powers of (j)/H, which would 
imply that the coupling to the spacetime metric (and all backreaction from it) switches 
off order-by-order in the limit <fi/H —* 0, and therefore that the field fluctuations are 
purely Gaussian in the de Sitter limit. Our result for M.± shows that this does not 
happen: the scalar fields and the metric remain coupled even in exact de Sitter space, 
and the fluctuations are therefore non-gaussian. Intuitively, this occurs because the 
size of the non-gaussianities in the 5(j) is controlled by the strength of the gravitational 
interaction, rather than (j)/H, even though powers of slow-roll appear once one has 
changed variable to the comoving curvature perturbation, (. We note, however, that £ 
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is not an appropriate choice of variable to discuss the limit <j)/H — * 0, since it becomes 
singular there. 

In the case of single-field inflation, we have applied Maldacena's consistency 
argument to the trispectrum. For the first time, this has employed the third-order term 
in the expansion of 5N in terms of initial field values. Only the linear and quadratic 
terms in the expansion have previously been required. In the case of the trispectrum, 
the consistency relationship arises as a next-order effect in the slow-roll parameter e, in 
contrast to that of the bispectrum. 

We have also derived expressions in the spatially flat gauge for the second-order 
metric perturbation sourced by multiple scalar fields. These expressions show explicitly 
that a metric perturbation continues to exist at this order even when the background 
is pure de Sitter space, corresponding to e — > 0. At first-order, the metric perturbation 
disappears in this limit. 
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